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Abstract. This paper is to explore the de Haas-van Alphen effect (dHvA) of 
graphene in the presence of an in-plane uniform electric held. Three major findings 
are yielded. First of all, the electric field is found to modulate the de Haas- van Alphen 
magnetization and magnetic susceptibility through the dimensionless parameter (/3 = 
As the parameter j3 increases, the values of magnetization and magnetic 
susceptibility increase to positive infinity or decrease to negative infinity at the exotic 
point j3 c = 1. Besides, the dHvA oscillation amplitude rises abruptly to infinity for zero 
temperature at /3 C = 1, but eventually collapses at a finite temperature thereby leading 
to the vanishing of de Haas- van Alphen effect. In addition, the magnetic susceptibility 
depends on the electric and magnetic fields, suggesting that the graphene should be a 
non-linear magnetic medium in the presence of the external field. These results, which 
are different from those obtained in the standard nonrelativistic 2D electron gas, are 
attributed to its anomalous Landau level spectrum in graphene. 
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1. Introduction 

Owing to the progress of experimental methods, graphene (or a graphite monolayer) is 
now gaining increasing interests in the field of physics of electronic systems with reduced 
dimensionality [1-4]. It is promising to be applied in nanoelectronics because of the 
exotic chiral features [5-9] in its electronic structure. In particular, such two-dimensional 
(2D) or quasi- two- dimensional systems have led to some of the most startling discoveries 
in condensed matter physics in the recent years. Moreover, these anomalous phenomena 
are found to be tied to the remarkable "relativistic-like" spectrum of electrons and holes 
in graphene, which makes graphene important and interesting in physics. One of them 
that have been experimentally testified is the abnormality of the 2D quantum Hall 
effect [10,11]. 

Another important physics effect is the de Haas-van Alphen (dHvA) oscillation of 
the graphene. It has been predicted in Ref. [12], which proposes that the magnetization 
oscillates periodically in a sawtooth pattern as a function of 1/B, in agreement with the 
old Peierls prediction [13]. A question of great interest arising here is what will happen 
if an additional electric field is applied in graphene. Indeed, the electric and magnetic 
fields effects on its magnetization and magnetic susceptibility are of vital significance 
to our understanding of the Dirac fermion behaviors. However, little theoretical or 
experimental research has been done on this issue yet. 

Motivated by the concerns mentioned above, the present study is to investigate 
the 2D dHvA effect of the graphene in the presence of the electric field. The paper 
is organized as follows. In section 2, a brief introduction is given to the 2D model 
for graphene. The energy eigenvalues and eigenstates, as well as the density of states 
(DOS) and the dHvA oscillation period A(l/B) are obtained analytically. Section 3 
describes some details of the magnetization and magnetic susceptibility study. In the 
meantime, the analytical expressions are derived for the magnetization and magnetic 
susceptibility. More specifically, the one regarding the condition of T = K is reported 
in section 3.1 while the one pertaining to the condition of nonzero temperature in section 
3.2. The section winds up with a discussion of the modulation of dHvA in graphene by 
the electric field. In the last section, the conclusions are presented. 



2. Energy eigenvalues and eigenstates 



In order to investigate the modulation of dHvA, we begin with the study of the 
energy eigenvalues and eigenstates belonging to the carriers in graphene. The charge 
carriers in graphene mimic relativistic particles with zero rest mass and have an 
effective 'speed of light' c* = v F « 1.0 x 10 6 ms" 1 , which is essentially governed 
by Dirac equation [10, 11]. So we start by considering the Dirac equation for such a 



2D gas of Dirac fermions in crossed electric 



E = (-E, 0, 0) , U = Eex 



and magnetic 



B = (0,0, B), A = (0,Bx,0) 



fields, where E is the electric field strength and B the 



The Modulation of de Haas-van Alphen Effect in Graphene by Electric Field 



3 



magnetic induction intensity. The single particle Hamiltonian is then given by 

H = v F a ■ n + leEx (1) 

in which a is the Pauli matrix, II is the canonical momentum, / is the 2x2 unit matrix. 
Following Landau and Lifshitz [14] the first-order equation of eigenvalue problem of H 
becomes the second-order equation 



e — eEx) — ( vpft '— eA ) + ehBv F a z + iehEvp6t x 



# = 0, 



(2) 



where e is the eigenvalue of H and other notations are standard. From Eq. (2), we 
obtain the energy spectra and eigenfunctions of the problem, 

in) ^2 \n\ heBv F (l - /3 2 ) 3/4 + Hv F /3k y 



£n,k y = sgn 



x y > 



with 



x + l 2 c k y — sgn(n) 




(1 - /3 2 ) 1 / 4 



(3) 
(4) 

(5) 



In Eq. (3), e and h = h/2n are electron charge and Planck's constant divided by 
27r, respectively. The integer n represents the Landau level index, k y = 2irl/L y 
(I = 0, ±1, ±2, • • •, ) is the quantum number corresponding to the translation symmetry 
along the y axis, L y stands for the size of the graphene in y direction. The electric 
field dependent dimensionless parameter (5 is defined by (5 = Ej (v F B) and obeys 
\/3\ < 1, where v F is the Fermi velocity. In Eq. (4), n (£) are the harmonic oscillator 
eigenfunctions. From Eq. (5) we observe that the centers of the x-dependent orbits are 
located at 



x = l c k y — sgn(n)- 



(6) 



;i _02)i/4 

where l c = ^Jh/eB is the magnetic length. The eigenvalues of H show that the exact 
energies are given by the sum of quantized harmonic-oscillator energies and the potential 
energy of a charged particle located at coordinate xq in potential field U(x). They agree 
with those in Ref. [15], the authors of which solved the problem by transforming the 
original system into a case with the null electric field, in terms of a Lorentz boost 
transformation. 

We then count the Landau states ^ n ,k y m the presence of the potential U(x) 
following the same argument employed in the absence of crossed electric field. Since 
k y = 2nl/Ly, the separation between adjacent allowed k y values is given by Sk y = 2n/L y . 
From Eq. (6) we may relate the possible range Ak y of k y to the physically accessible 
range Ax of x : 

Ax = l\Ak y . (7) 
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Since Ax = L x , in order for the Landau states to be centered within the strip 
< x < L x we must have allowed the range of k y values given by Ak y = Axq/1^.. 
The number of Landau states ^ n ,k y (x, y) per unit area for each quantum number n is: 

Dn = 9sAk y = (8) 

which is independent of n and the degeneracy yields g s = 4, accounting for spin 
degeneracy and sublattice degeneracy in graphene. 

We consider a system of iV electrons within an area S moving in the potential 
U (x) and the magnetic field B. Let the system remain at K and accordingly the free 
energy reduces to the total energy. The full occupation of Landau levels obeys D n S = 
N/ (n F + 1) with Fermi quantum number n F . The total energy E = ^2 nky S n ,k v (B) 
will give a discontinuous derivative M = —OE/dB at the field values B n where M is 
the magnetization. From D n S = N/ (np + 1) these discontinuities in the magnetization 
occur at reciprocal fields l/B n , so that the period of magnetization oscillation is given 
by: 

1\ 1 1 2e 

(9) 



B) B n+l B n nhN 

where B n+i and B n are the magnetic induction intensity corresponding to two 
neighboring levels, which cross the Fermi level in succession and N = N/ S is the sheet 
concentration. Eq. (9) means that the discontinuous zero-temperature oscillations are 
periodic in 1/B. It is just compatible with the results regarding the null electric field 
obtained by Sharapov et al [12]. 



3. Magnetization and magnetic susceptibility 

3.1. Zero temperature 

Then we move on to investigate the magnetization of electrons in graphene in the 
presence of crossed uniform electric and magnetic fields at T = K. For simplicity, 
we ignore spin-orbit coupling of electrons in the present work. The magnetization 
reads [16, 17] 

M — — (dE/dB) N (10) 

where E is the total energy and iV denotes the total number of electrons in graphene. 
To have units in Tesla, we symbolize B = n P oH, /x p o as being the magnetic permeability 
of free space, H stands for the magnetic field intensity. 
The total energy is given by: 

[n%/2ehv F B\ 

E = Y1 En ^y = Dn ( £n ~ V) + V N o 

n,k y n=0 

+ J2 ^v F pi,(-i F <i y <i F ) (ii) 
i=~i F Ly 
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in which the last term is the additional energy induced by the electric field corresponding 
to the npth level partly occupied by electrons and we refer to it as E a &&. Here e n is 
defined clS En 



y/2nheBvp (1 — /3 2 ) 3 ^ 4 . The chemical potential ji 



Mi-/? 2 ) 3/4 



is 

derived analytically where fio refers to the zero-temperature chemical potential (equal 
to the Fermi energy) in the absence of electric and magnetic fields as expressed by 
/io = hvFy/Noir. Using the formula for the generalized zeta function [18], ( (z, v + k) = 
C v) - Z k ~=o ( m + V Y Z and C (-1/2, 0) = C (-1/2) = - (1/4tt) C (3/2), one can write 
the first and the second terms in Eq. (11) as the summation of the regular term, 



E, 



C (3/2) v F 



eBf 2 (l-/3 2 )^ + 



2\3/4 



reg 



7T 



2h 



3n {Hvf)' 



2\3/4 



(12) 



and the oscillating term 



E SC 



2V2hv r 



(eBf 2 (1 - (3 2 ) 



irh 



2\3/4 



C[- 5 ,i + 



^0 



2ehv 2 F B 



- (i - f) 



2\3/4 



2/x 



2fi eB 
3tt (hv F ) 2 



mod 



t4 



2ehv F B 



(13) 

In this expression, [^/2ehv F B] stands for the integer part of $/2ehv F B and the mod 
\p^/2ehv F B\ is the fractional part of ^/2ehv 2 F B. 

Making use of Y function r (n + a) = J °° dss n+a ~ 1 e~ s and Bernoulli polynomials 



£ 

n=2 
f oo 



r (n + a) B n (y) 



x' 



TV. 



dss a ~ x e- s 



sxe 



sxy 



- 1 



1 - sxB x (y) 



(14) 



where the explicit expressions for the Bernoulli polynomials B , Bi, B 2 are 

B (y) = 1, B 1 (y) = y- 1/2, B 2 (y) = y 2 - y - 1/6. (15) 

The Bernoulli polynomials depend on mod[y] in the following equations, i.e. 
B n (mod [y]) . For brevity, we write it as B n (y). 
Using the formula [17] 



oo x v-l e -»x dx 



1 - e-P* 



(3 V 



i-F (v)((v,^),(Rev>0,Rev>0) 



we have 



2 (eBv F ) 

7T 3 / 2 /i 



2 , 3/4 A T (n + 1/2) B n+2 (w/2) (2ehv 2 F B 



?1=0 



(n + 2)\ 



t4 



(16) 

in which w = //jj/ (eh~v F B). Regarding small fields, eBhv F <C //q, we can apply the 
following asymptotic expansions for Ji = / °° dte~ tp / [y/nt (t 2 + 1)] : 

-l) w r(n + l/2) 



■AC?) 



n=0 



(17) 
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and the Bernoulli polynomials B n periodically continue beyond the interval [0,1]: 

2n\ 



B n = 



(27T) 



k=i 



(18) 



n 



> 1,0 < x < l;n = 1,0 < x < 1. 
It is easy to get the following form: 



E = 



{eB) 3/2 v F 
\ffm 



(i - p 2 ) 



k=l 



1 



[irk) 



3/2 



J\ (irkw) cos (nkw) . 



(19) 



For \J 'eBhvp <C /i , keeping the leading term in asymptotic expansions for J\ (p), we 
finally obtain from Eq. (19) 



E — 

^osc — 



{eBv F f 

7T/i 



cos (irkw) 



(20) 



Hence the total energy E can be expressed as a sum of regular, oscillating and the 
additional energy terms, 

E = E reg + E osc + E add . (21) 

According to the results reported above, we get the corresponding de Haas- van 
Alphen magnetization, 

3eVeBh<; (3/2) v F /ifj/3 2 



M reg = 



M, 



2V2n 2 h (1 - /3 2 ) 1/4 7T (hv F ) 2 B (1 - /3 2 ) 1/4 ' 
irhv F {3 



add 



BL, 



-A 



M 5 



where 



/4 



(^f) 



7T 



/*0 



2ehv F Bj 



t4 



2ehv 2 F B 



+ 1 



2eB 



(22) 
(23) 

(24) 



This expression involves a dependence on the integer part [fi 2 ) /2ehv F B] and 



M n 



A 



Ml 



oo 

E 

fe=i 



7r/c 



sin 



/ nk^l \ A ^ -1 f f 

r ^ + A M2 > 1 COS 



nkji 



ehvpB 



(25) 



where A M1 = 



3/4 



7T/1 



and A 



M2 



(ev F ) 2 B(4-/3 2 ) 
' 2tt M o(1-/3 2 ) 1/4 



Figure 1 shows that the dHvA oscillation on the magnetization M modulated 
by electric field in graphene. As shown in Fig. 1(a), the magnetization M oscillates 
periodically in 1/B with the period of A (1/5) = 2ehv 2 F /n 2 ) . Three oscillation curves 
correspond to the applied electric field strength E 1 — 5 V/m, E 2 = 10 V/m and E 3 = 15 
V/m, respectively. From these curves, one can see that the oscillation amplitude {OA) of 
magnetization is proportional to the electric field E. Also, the OA ascends significantly 
with increasing 1/B, but remains unchanged for the null electric field in Ref. [12]. 
Thus, it is demonstrated the electric field effect on the OA of magnetization. However 
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FIG. 1 . (Shengli Zhang et al) 
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Figure 1. (a) Magnetization M is plotted as a function of reciprocal magnetic field 
1/B for a given chemical potential /xq = 1-6 x 1CP 19 J and T — K. Three oscillation 
curves correspond to E\ = 5 V/ m, E% = 10 V/m and E 3 = 15 V/m, respectively. (6) M 
plotted vs electric field E for /xq = 1.6 x 10 -19 J and T — OK. Three curves correspond 
to B\ = 2 T, i?2 = 3.63 T and B 3 = 5 T, respectively, (c) The magnetization, M, as 
a function of 1/B at E = 0.001 V/m, ^ = 1.6 x 10~ 19 J and T = K. 



unexpected it may be, according to Eq. (25) the OA increases abruptly to infinity 
at 0c = 1 thereby leading to the vanishing of dHvA oscillation. It is interesting that 
the abnormal phenomena will die out if the electric filed vanishes. Thus, it can be also 
attributed to the electric field effect. We hope that the new findings will be verified when 
magnetization experiments under in-plane electric field are carried out in graphene. The 
predicted effect will hopefully also help interpretation of magnetization in experiment. 

Meanwhile, we find that these peculiar features are absent in standard quantum 
2D electron gas systems [19-21]. Hence, the possible reason for the effect might be that 
it is determined by the "relativistic" character of carriers in graphene, unlike the usual 
sample, which can be traced to the exotic structure in graphene. 

Figure 1(b) demonstrates that the magnetization M is a function of the electric field 
E. Three curves correspond to B 1 — 2 T, B 2 = 3.63 T and B 3 = 5 T, respectively. They 
exhibit the magnetization varies approximately linearly with increasing E within the 
given values of parameters. For B\ = 2 T, the magnetization increases with increasing 
E. But for B 3 = 5 T, as E increases the magnetization decreases. Especially, there is 
a special behavior of the magnetization at some magnetic fields (e.g. B 2 = 3.63 T). 
In this case, the magnetization satisfies M ~ 0, accounting for the disappearance of 
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magnetization in graphene. Note that, the magnetization become infinite when the 
variation of E obeys ft c — 1 according to Eqs. (22) and (25). For example, the dashed 
line (B\ = 2 T) rises abruptly to infinity at E = 2 x 10 6 V/m, strikingly different from 
the non-relativistic results. 

In Fig. 1(c), we present the dHvA effect for a wide range of magnetic fields starting 
from 10 to 400 T, in order to examine what will happen if the magnetic field tends to 
be infinite. As a result, we observe that on this scale, the electric field effect could be 
negligible corresponding to the case ft — > 0. In other words, for the case of ft — )■ 0, the 
magnetization M reduces to the result for the null electric field. 

The same discussion above fits for the magnetic susceptibility x- We can obtain 
the expression of the dHvA magnetic susceptibility in terms of x — dM/dH, 

/ig/3 2 /V) (6 - 5ft 2 ) 



and 



Xreg 



Xadd 



Xosc 



3e//poVefa(3/2) v f / _ ^ 

32)5/4 I P ) 



W2BTT 2 h(l 

_ 2/i p0 /3^Q 
L y B 2 (hv F ) 



2tt (hv F By (1 - ft 2 ) 



2^5/4' 



3 ' 



(26) 



(27) 



(A xl + A x2 )J2(-±) k c°s 
+ {A x3 + A 



k=l 



ehvpB 



k=i 



7ik 



sin 



ehvpB 



(28) 



where 



2N3/4 



A 



X3 



7T (HvfB) 

UpoefiQ (4 - ft 2 



1/4' ^X4 



(ev F ) 2 fi p0 (8 -10 ft 2 -ft 4 
4n^ (nk) 2 (l-ft 2 f 4 
3/i p0 efi ft 2 



(29) 

2nhB{l- ft 2 )' 1 *' A ' 2nhB (1 - ft 2 ) 1/4 
We can see that the magnetic susceptibility is related to both the electric and magnetic 
fields. It follows that, unlike the usual samples, graphene may be a non-linear magnetic 
medium. 

Figure 2(a) shows that the magnetic susceptibility x oscillates periodically as a 
function of 1/5 and the period follows Eq. (9). Three oscillation curves correspond to 
E 1 = V/m, E 2 = 2 V/m and E 3 = 4 V/m, respectively. For the case (Ei), the reader 
may observe that the magnetic susceptibility x swings between negative and positive 
values, thus the curve shows a totally orbital diamagnetic to paramagnetic transition. 
As for the finite electric fields, one can see that the magnetic susceptibility decreases in 
company with the periodic oscillation while 1/B rising and the OA augments as 1/B 
increases. Furthermore, from Eq. (28) it can be seen that the OA increases to infinity at 
ft c — 1, leading to the vanishing of dHvA effect on magnetic susceptibility. In general, 
the magnetic susceptibility x is a constant in usual electron gas, but in graphene it 
exhibits the dependence on the external field. 
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FIG. 2. (Shengli Zhang et al) 
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Figure 2. (a) Magnetic susceptibility \ depends on reciprocal field 1/B for a given 
chemical potential /j,q = 1.6 x 10~ 19 J and T = K. Three oscillation curves of 
magnetic susceptibility correspond to E\ = V/m, E% = 2 V/m and £3=4 V/m, 
respectively, (b) \ as a function of electric field -E for fiQ — 1.6 x 10~ 19 J and T = OK. 
Three curves correspond to B\ = 2 T, B2 = 3 T and B3 = 15 T, respectively. 

Figure 2(b) depicts the magnetic susceptibility x with respect to the electric field 
E. Three curves correspond to B\ = 2 T, B2 = 3 T and — 15 T, respectively. As 
is shown by the three curves, the magnetic susceptibility varies approximately linearly 
with increasing E. For B = 2 T, or 3 T, the magnetic susceptibility decreases as 
E increases and yields x < 0, indicating the existence of Landau diamagnetism in 
graphene, the origin of which can be traced to the quantized Landau level. In the case 
of B = 15 T or larger, the magnetic susceptibility \ ~ despite the increase of E, 
suggesting the disappearance of the diamagnetism in graphene. That is, there is no 
increase in the magnetic susceptibility with increasing E. Moreover, the dashed line 
{Bi = 2 T) decreases to negative infinity at the exotic point E = 2 x 10 6 V/m (i.e. 
p c = 1) as illustrated by Eqs. (26) and (28). 

3.2. Finite temperature 

We now consider the temperature effect on the oscillations of magnetization and 
magnetic susceptibility. As documented in [12], the thermodynamic potential of 
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/oo 
duP T (u-n)E (to) , 
-oo 

with energy variable u> and /i. And P T (z) is the distribution function as 

dnp (z) 



P T (z) 



dz 



1 

4/csTcosh 2 



2k B T 



Using Eq. (30), the thermodynamic potential Q can be divided as follows: 

^ (-^"'5 AV ^reg ~t~ ^add ~t~ ^osc- 

At the low temperatures, neglecting 0(k B T), we can obtain 



C (3/2) v F 



(eBf 2 (l-^Y / ' l + 



2\3/4 



reg 



2/xg 



37T (hVp)' 



(1 - 



2\3/4 



7T 



^add = —hv F /3A nT R T (h,(i), 
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(30) 
(31) 
(32) 

(33) 
(34) 



where 



An, 



t4 



2eBfi 2 
7r 2 [hvpY hinhvp) 



t4 



2ehv 2 F B 



+ 1) + 



tA 



2ehv 2 F B 



x 



tA 



2ehv 2 F B A 



(2eB 



(35) 



This expression involves a dependence on the integer part [^/2ehv F B]. And here we 
introduced the temperature reduction factor 

27v 2 kj2 k B T/ (ehv 2 F B*) 



R T (k, 11) 



sinh 2 ^ k R f 



(36) 



in which B* = B (1 — (3 2 ) 3 ^ 4 . Eq. (36) means that the temperature reduction factor Rt 
depends not only on the temperature, but also on the electric and magnetic fields. We 
return now to the oscillating part of thermodynamic potential. Substituting Eq. (19) 
into Eq. (30), it is convenient to get the expression of fl osc as 



x Im [e 



-17T/4 



(^) 3/2 

dte -i(nkw)t 



de 



y/i(t + l)J-oo Ak B T cosh 2 §g£ 



-mke 2 (t + 1) 



^Wfj 1 (37) 

in which we used that the function J±(p, r) and J 2 (p, f) which can be represented as -Im 
and Re parts of the same function 

dte- pt ~ r/t 



\fnJi (p, r) = — Im 
\fwJ 2 (p, r) — — Re 



y/i(t + i)' 
00 dte- pt - r/t 
Vt(t + i)' 



(38) 
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and rotated the integration contour to the imaginary axis. Finally we get 



«~ (T, „) = (1 - £ ^^Rt (*, M) , (39) 

7T/X ^ (7Tfc) 

with iu = nl/ (ehvpB). Clearly, since R T (k,fx) -> 1 for T — >■ 0, Eq. (39) reduces 
to the oscillating energy for zero temperature. When the temperature T / 0, the 
magnetization can be obtained as follows: 

T 3eVemC (3/2) v F _ ggg 

reff 2^2^(1-^)1/4 7r(ftujil) 2 s(1 _ ^ 2) V4' ^ ^ 

T 7rhv F f3R T 7rhv F f3 (2 + /3 2 ) 



and 



6/i ° Cl g2\ 3/4 sin (7rfcw) _ _{evpfB__ 
nh [ ' ' ^ nk T 27r/i (l-/3 2 ) 1/4 



fc=i 

x [(4 - /3 2 ) i? r + (2 + /3 2 ) (i$ - ifc)] £ (42) 

where R F (k,fi) = R\ (k, /i) cosh (2ir 2 kfj,okBT /ehvpB*). Apparently since Rr(k,fi) — > 
1 and Rtp (k, fi) — > 1 for T — > 0, Eq. (42) reduces to the magnetization in Eq. (25) for 
zero temperature. 

As is illustrated in Figure 3(a), finite T causes a reduction of the magnetization 
amplitude as opposed to the case of T = K. Meanwhile, it shows the magnetization 
M versus reciprocal field 1/B for three different temperatures with a given E. It can 
be seen that as the value of 1/B increases, the OA of magnetization decreases and 
eventually collapses when j3 — > 1. Nevertheless, regarding the zero temperature, there 
is no such a collapse for the OA of magnetization as shown in Fig. 1(a). Accordingly, 
we attribute them to the finite temperature effect. 

Figure 3(b) depicts the magnetization M as a function of E for three different 
temperatures with a given B. Three curves correspond to 7\ = 20 K, T 2 = 30 K 
and T3 = 100 K, respectively. It has been shown that the magnetization decreases 
approximately linearly with increasing E except in the case of T 3 = 100 K. Finally, 
the magnetization decreases to the negative infinity at (3 C = 1 following Eqs. (40), (41) 
and (42). At the temperature T 3 or even a higher one, the magnetization M ~ rather 
than decrease with increasing E. Fig. 3(c) gives the magnetization M as function of 
the temperature T. Increasing T, the magnetization M increases for Bi, B 3 , B 5 and 
decreases for B 2 , B±, B 6 , but they all end up approaching different constants. 

We get the dHvA magnetic susceptibility from the magnetization: 

t 3efi p oVe~hq{3/2) v F , _ 2 . _ /ig/3 2 /i p0 (6 - 5(3 2 ) 



AV2B7r 2 h (1 - /3 2 ) 5/4 V ' 2tt {hv F Bf (1 - (3 2 f 4 ' 
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FIG. 3. (Shengli Zhang et al) 
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Figure 3. (a) Magnetization M is plotted as a function of reciprocal magnetic field 
1/B for a given chemical potential /io = 1.5 x 10~ 19 J and E = 2 V/m. Three curves 
correspond to temperature T\ = 90 K, T 2 = 100 K and T 3 = 150 K, respectively, (b) 
M plotted to electric field E with B = 2 T and fi = 1.5 x 10~ 19 J. Three curves 
correspond to T x = 20 K, T 2 = 30 K and T 3 = 100 K, respectively, (c) M is plotted 
as a function of temperature T for /io = 1.5 x 10~ 19 J. 
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/(.add 
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and 



T 

Xosc 



v4 x i T cos (nkw) + A X 2 T sin (nkw) , (45) 
fc=i fc=i 
in Eq. (44), A Mo and Aq t are denned by Eqs. (24) and (35), respectively. We also 
define, 

. (2/3 4 - 10/3 2 - 4) / 2 + /3 2 \ 2 

Ax ° T ~ [25(l-/3 2 )] 2 [Rt ~ Rt]+ \2B{1-^)) [2Rt 

(46) 
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xlT ~ ^ — m2 ^ .2,. ^TsTi — v t - ^tJ 

7r (Kv F B) 2ttij, (irk) (1 - f3 2 ^ ' 



The Modulation of de Haas-van Alphen Effect in Graphene by Electric Field 



13 



FIG. 4. (Shengli Zhang et al) 





Figure 4. (a) Magnetic susceptibility x to reciprocal field 1/B for a given chemical 
potential /io — 1.5 x 10 -19 J and E = 0.01 V/m. Three oscillation curves of magnetic 
susceptibility correspond to T\ = 80 K, T 2 = 100 K and T3 = 150 K, respectively. 
(b) x as a function of electric field E with B = 2 T and fi = 1.5 x 10~ 19 J. Three 
oscillation curves of magnetic susceptibility correspond to T\ = 20 K, T 2 = 30 K and 
T3 = 100 K, respectively, (c) \ as a function of temperature T for /iq = 1.5 x 10~ 19 J. 



_ R (et; F ) 2 /i p o(8 - 10/3 2 - /3 4 ) _ ^ x o T /x p0 (ev F B) 2 , _ 2 , 3 /4 
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(47) 

A x2t = -^^±4,^. (48) 
XT n 2 hkB(l -/3 2 ) 1/4 

Figure 4(a) shows the magnetic susceptibility \ oscillates periodically as a function 
of 1/5, for three characteristic temperatures, namely, 80, 100, and 150 K. Similar to the 
magnetization as shown in Fig. 3(a), it also exhibits the dependence on temperature. 
It can be seen clearly that finite T causes a reduction of the oscillation amplitude and 
as the value of 1/B increases, the OA of magnetic susceptibility finally decays to zero. 

Figure 4(b) shows the magnetic susceptibility \ plotted as a function of E with a 
given B. Three curves correspond to T\ = 20 K, T2 = 30 K and T3 = 100 K, respectively. 
For a finite temperature such as T\ or T 2 , the magnetic susceptibility decreases 
approximately linearly with increasing E and exhibits the Landau diamagnetism in 
graphene supported by x < 0. In contrast, at the temperature T 3 , the magnetic 
susceptibility x — 0, standing for the disappearance of the diamagnetism in graphene. 
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In addition, there is no increase in the magnetic susceptibility with increasing E. It 
warrants great attention that from Eqs. (43), (44) and (45), it could be inferred that 
the magnetic susceptibility eventually decreases to the negative infinity at (3 C — 1. Fig. 
4(c) directly gives the magnetic susceptibility x with respect to the temperature T. 
More specifically, it shows that the magnetic susceptibility x increases with increasing 
T and finally approaches zero under different electric field Es and the magnetic field 
Bs. 

4. Conclusions 

In summary, this paper reports on a theoretical study on the modulation of de Haas- van 
Alphen effect in graphene by electric field. Three major findings emerge from the study. 
First of all, we find that both magnetization and magnetic susceptibility are modulated 
by the electric field. At the zero or finite temperature, both magnetization and magnetic 
susceptibility are predicted to oscillate periodically as a function of reciprocal field 1/B. 
The dHvA oscillation period A(l/B) is derived analytically. It is also discovered that 
as the parameter j3 increases, the values of magnetization and magnetic susceptibility 
finally increase to positive infinity or decrease to negative infinity at the exotic point 
(3 C — 1. Besides, the analytical results indicate that the dHvA oscillation amplitude 
increases abruptly to infinity for zero temperature at (3 C — 1, but eventually collapses 
at a finite temperature directly leading to the vanishing of the de Haas-van Alphen 
effect. The "vanishing" is accounted for the anomalous effect of the electric field on 
Landau level, which arises from the instability of a relativistic quantum field vacuum. 
In addition, it is established that the magnetic susceptibility depends on the electric 
and magnetic fields, which suggests the graphene should be a non-linear magnetic 
medium. These phenomena, not available in the standard 2D electron gas, are deemed 
as the consequence of the relativistic type spectrum of low energy electrons and holes 
in graphene. 
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